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ACTIVE GALACTIC NUCLEUS PAIRS FROM THE SLOAN DIGITAL SKY SURVEY. 
EVIDENCE FOR TIDALLY ENHANCED STAR FORMATION AND BLACK HOLE ACCRETION 
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ABSTRACT 

Active galactic nuclei (AGNs) are occasionally seen in pairs, suggesting that tidal encounters are responsible 
for the accretion of material by both central supermassive black holes (BHs). In Paper I of this series, we 
selected a sample of AGN pairs with projected separations r p < 100 h^,' kpc and velocity offsets < 600 km 
s" 1 from the SDSS DR7 and quantified their frequency. In this paper, we address the BH-accretion and recent 
star-formation properties in their host galaxies. AGN pairs experience stronger BH accretion, as measured by 
their [Om] A5007 luminosities (corrected for contribution from star formation) and Eddington ratios, than do 
control samples of single AGNs matched in redshift and host-galaxy stellar mass. Their host galaxies have 
stronger post-starburst activity and younger mean stellar ages, as indicated by stronger HS absorption and 
smaller 4000 A break in their spectra. The BH accretion and recent star formation in the host galaxies both 
increase with decreasing projected separation in AGN pairs, for r p < 10-30 h^ kpc. The intensity of BH 
accretion, the post-starburst strength, and the mean stellar ages are correlated between the two AGNs in a pair. 
The luminosities and Eddington ratios of AGN pairs are correlated with recent star formation in their host 
galaxies, with a scaling relation consistent with that observed in single AGNs. Our results suggest that galaxy 
tidal interactions enhance both BH accretion and host-galaxy star formation in AGN pairs, even though the 
majority of low redshift AGNs is not coincident with on-going interactions. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei - galaxies: 
stellar content 



1. INTRODUCTION 

It has long been recognized that galaxy mergers and close 
interactions pla y a vital role i n shaping the morphology of 
galaxies (e.g., IZwickvl [T 956: Vorontsov-Velvaminov 1 19591 
Toomre & Toomrd 119721: iKhachikianl 11987b ISanders et al.1 
1988t iBarnes & Hernquistl 1 1992b and in particular the 



cores of ellipticals (e.g., Faber et all [l997t iGrahaml 120041 



ellipticals 

Merrittetail l2004t lLauer et al.1 120051 iKormendv & Bended 
20011 7 While mergers only represent brief episodes in the life- 
time of a galaxy (e.gjLotz et alJl2008HConselice et al]l2009l : 
Darg et al. 2010bt lLotz et al.l201 ll) . there is circumstantial ev- 
idence that they can trigger substantial star formation both in 
the nuclea r region and throughout the galaxy (e.g., see the 
review by Struck 2006). The enhancement in star forma- 
tion predicted by numerical simulations depends on how that 
star formation and its associated feedback are modeled (e.g. , 
Mihos iraTJ[T99l iK auffmann et all [1991 ICole et alj \l994 . 
Sprin gel & Hern quist 2003; Springel et al. 2005b). Therefore 



observations of star formation in interacting galaxies can be 
used to constrain these models, and help advance our general 
understanding of star formation, and in partic ular, modes that 
are exclusive to strong t idal encounters (e.g., Elmegree n et alJ 
fT99llDaddTet alJ2010l) which may be more commo n in early 
galaxy formation (e.g jForster Schreiber et al.1l2009b . 

Galaxy mergers and close interactions are also thought to 
drive significant ac cretion in their central supermassive black 
holes (BHs; e.g.. iHernquisil 119891; iMoore et al.1 1 19961) and 
are thus responsible for at least some Active Galactic Nu- 
clei (AGNs). The enhanced nuclear activity may be pri- 
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marily driven by tidal torques, or secondarily fueled by stel- 
lar mass loss due to elevated levels of star formation (e.g., 
Norman & Scovillell988T:ICiotti & Ostrikerl2007UCiotti et all 
20101) . Some numerical models predict that merger-induced 
AGNs heat the gas surrounding th em, preventing further 
gas cooling and star formation (e.g ., Springel et al. 2005b a; 
iDi Matteo et alj|2005t iHopkins et alj|2006l) . Observationallv. 
however, the causal link between galaxy interactions and 
AGNs has been controversial. The details of AGN feed- 
ing, and in particular, the role of mergers, ma y be quit e 
different at high and low accretion rates (e.g., iHol 2008). 
While the small-scale (~ 0.1-1 Mpc) quasar-quasar two- 
point correlation function suggests a c lustering excess over 
the large-scale (> 1 Mpc) extrapolation (IHennawi et all2006fc 
iMvers et al]l2007t IHennawi et alll20Tot Ishen et al]l201TbTu t 
is unclear whether this is due to interactions of the quasar 
pairs' host galaxies (e.g.. |Diorgovskilll99lt iKochanek et al. 
119991: iMortlock et al.l 1 19991) . or is rather a manifest of the 
small-sc ale clustering strength of their host dark matter ha- 
los (e.g., IHopkins et alJl2008h . Some observations found an 
excess of close neighbors in AGN host galaxies or a higher 
fraction of AGNs in interacting than in isolated galaxies (e.g. , 
Petrosianlll982l: iHutchings & Campbell 1983b [Paharil 1 1984 : 
Kennicutt & Keell 1984l:lKeel et alJI 1985!:lBahcall et al] 1997 : 
Serber et al.l ) 2006 t iKoss et all l2010t ISilverman et al.l boi 1 



Ellison et al.ll201 1|) wherea s others detected no significant dif- 
ference (e.glDaharill 19851: ISchmitdl20qil iMiller et al.l | 2003 
i Grogin et al] 120051; IWaskett et al 1 120051; Ellison et all 12008 
Ili et al]|2008btlDarg et al.ll2010bl) . 

Observations of galaxy pairs directly probe the effects in- 
duced by galaxy tidal interactions, at least at the early stages 
when the two galaxies can be identified as separate. Studies 
of the properties of galaxy pairs as a function of pair sep- 
aration, and in comparison to single galaxies provide clues 
to how tidal effects evolve as mergers progress. Simulations 
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FIG. 1 . — Redshift and stellar mass distributions of AGN pairs (left panel; blue) and the subset with tidal features (right panel; blue) compared to their control 
samples (magenta in both panels). The AGN pair and control samples are matched in redshift and stellar mass (with KS probabilities larger than 50%). The 
median value of each distribution is indicated with a vertical line. 



suggest that merger-induced star format ion and AGN activity 
peak after the galaxies have coalesced ( Sp ringel et i 



jfN activity 
iD|2005b; 



iDi Matteo et al.l 120051) . There is accumulating evidence for 
tidally-enhanced activity starting at the early merger stages 
when the galaxies are still separate and that the enhancement 
becomes more prominent as a merger progresses. Emission- 
line galaxies in pairs show stronger star format ion activity. 



which is larger for smaller projected separ ations (Barto n et al 
2000t Lambas et al. 2003t lAlonso et ai.ll200H iNikolic et al 



2004t Ellison et all 12008). Galax y-AGN pairs with smaller 
separations show stronger AGN ([Rogers et al.l 120 09), while 
galaxy p airs with smaller separations host a high er fraction 
of AGN dWoods & Gelled l2007t iDarg et al.ll2010ah . If tidal 
effects impact both interacting components (e.g., by promot- 
ing gas infall and inducing gravitational instabilities) and 
the induced activity is synchronous, there should be observ- 
able correlations between th em, unless additio nal factors reg- 
ulate such trends. Indeed, Holmberg (119581) remarked on 
what is now called the "Holmberg effect", w hereby the col- 
ors of galaxies in pairs are correlated (e.g., lMadorelll986[ 
iKennicutt et alJl987klLaurikainen & Moles! 19891) which sug- 
gests the presence of tidally induced star formation in both of 
the interacting galaxies. 

We have selected a sample of 1286 pairs of AGNs with 
projected separations r p < 100 h™ kpc and line-of-sight ve- 
locity offset s Av < 600 km s" 1 fr om the Seventh Data Re- 
lease (DR7; lAbazaiian etal]|2009h of the Sloan Digital Sky 
Survey (SPSS: lYork et al l (2000T) and examined their fre- 
quency among t he parent sample of optically selected AGNs 
dLiu et al .11201 ll Paper I). The sample is defined in detail in 
Paper I. In the present paper, we study their BH-accretion and 
recent star-formation properties in the host galaxies to con- 
strain the effects of galaxy tidal interactions. AGN pairs char- 
acterize a special population of galaxy pairs, in which the cen- 
tral SMBHs of both galaxies are actively accreting at the same 
time. A statistical sample of AGN pairs enables us to exam- 
ine whether there is a correlation in the accretion power of the 
central SMBHs between the interacting components, akin to 
what has been observed for star formation in galaxy pairs. 

In ^2] we describe the construction of our AGN pair sam- 
ple and the control sample of single AGNs. In §|3]we present 
the recent star formation and BH accretion properties of AGN 
pairs, and we examine their dependence on th e pro jected pair 
separation and recessional velocity offset in ^3.11 the depen- 
dence on host-galaxy properties in jj3.21 the correlation be- 
tween the two components in each AGN pair in ^3.31 and the 
cor relat ion between recent star formation and AGN activity 
in jj3.4l We discuss the implications of our results in ^4] and 
conclude in §|5] Throughout we assume a cosmology with 



ft,,, = 0.3, Q A = 0.7, and H = 70 h 10 km s" 1 Mpc 



-i 



2. DATA: AGN PAIRS, AGN PAIRS WITH TIDAL FEATURES, AND 
CONTROL SAMPLES 

We briefly describe the basic properties and physical mea- 
surements of our AGN pairs. Details of our sample selection 
are presented in Paper I. We have selected AGN pairs with 
r p < 100 hjQ kpc and Av < 600 km s" 1 from a parent sample 
of 138,070 AGNs spectroscopically identified based on their 
optical diagnostic emission-line ratios 4 and/or widths. 98% 
of the sample is co ntained in the SDSS main galaxy catalog 
(IStrauss et a l. 2002), and for the present a nalysis we have ex- 
cluded broad-line objects (lHao et al .120051) to avoid AGN con- 
tinuum contamination on host-galaxy measurements. We re- 
fer to this sample as the "pair" sample. As discussed in Paper 
I, 256 of the 1286 AGN pairs show unambiguous morpholog- 
ical tidal features tails in the SDSS images, such as bridges, 
and/or rings, indicative of strong tidal interactions, and we re- 
fer to this subset as the "tidal" sample. We consider the tidal 
sample to be "cleaner" than the pair sample for studies of in- 
teractions, as the pair sample may include closely separated 
AGN pairs that are not yet tidally interacting. On the other 
hand, as discussed in Paper I, the tidal sample is incomplete 
and subject to biases due to the surface-brightness and image- 
resolution limits in our ability to recognize tidal features. We 
thus analyze both samples and compare results to quantify the 
range of possible values. As discussed in Paper I, due to the 
finite size of the SDSS fibers, galaxy pairs with separations 
smaller than 55" will not both be obser ved unless they fall 
in the overlap regio ns on adjacent plates (IStrauss et al.ll2002T: 
iBlanton et al.l 120031) . We correct for this spectroscopic in- 
completeness by supplementing the observed sample of AGN 
pairs with (C- l)N p repeated systems randomly drawn from 
the AGN pairs with separations smaller than 55", where N p is 
the observed number of pairs and C w 3 .3 is the correction fac- 
tor (see Paper I for details). The small separation AGN pairs 
fall in the overlap regions on adjacent plates so that they both 
got spectroscopic observations. There should be no biases in- 
troduced here because the spectroscopically observed small 
separation pairs are randomly distributed on the sky and are 
therefore representative of the parent small separation popu- 
lation. 

We adopt redshifts z and stellar velocity disper- 



4 For narrow emission line objects, we adopt the Kauffmann et al. (2003a) 
empirical criterion for AGN classification as our baseline values (i.e., includ- 
ing AGN-H ii composites to avoid selection biases against galaxies with more 
starburst components) , but we have also t ested that our results do not vary 
qualitatively using the Kewley et al. 1 2001) theoretical "starburst limit" crite- 
rion (i.e., excluding AGN-H n composites to focus on AGN dominated sys- 
tems). 
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FIG. 2. — Stellar mass M t and r-band concentration index Rgtir/Rso.r as functions of interaction parameters r p and Av for the pair (left block) and the tidal 
(right block) samples. Here and throughout, we plot both galaxies in each AGN pair. Cyan dots denote individual objects, whereas blue filled circles represent 
median values at each bin of r p or Av. Error bars indicate uncertainties of the median values. The histograms plotted in thick curves in cyan are normalized 
distributions of AGN pairs corrected for SDSS spectroscopic incompleteness due to fiber collisions (Paper I). The thin dotted curves in magenta are for the 
control AGNs matched in redshift and in stellar mass (Figure [T}. Vertical and horizontal lines denote median values of the distributions. Here and throughout, 
we do not plot the median values below Av = 30 km/s [i.e., log(Av/km s -1 ) < 1.5] which are insignificant because they are dominated by redshift measurement 
uncertainties. 

sions a* of each galaxy from the specBS pipeline 
dAdelman-McCarthv et al.ll2008t; lAihara et aI1l201 It see also 
discussion in iBlanton et al.l 120051) and excluded extreme 
values of cr* (cr* < 30 or > 500 km s" 1 ) and those with 
large uncertainties (S/N< 3) or negative errors from bad 
fits. Additional spectral and photometric properties such 
as emission-line fluxes, continuum spectral indices, stellar 
masses, and half- light radii are take n from the MPA-JHU 
data product 5 (see lAihara et al.ll201 ll for a description of the 
catalog). The emission-line measurements a re from Gaussian 
fits t o continuum-subtrac ted spectra dBrinchman n et al.l 
2004; Tremontietafll2004l The emission-line fluxes have 
been normalized to SDSS r-band photometric fiber mag- 
nitudes and have be en correcte d for Galactic foreground 
extinction fol l owing lO'Donnelll (11994 using the map of 
Schle geTet al.l d!998f) . Continuum spectral indices such as 
H<5a and D„(4000) are measured from data spectra after 
subtracting all 3-cr emi ssion lines dKauffmann et alJ l2003bb 
iBrinchmann et al.l [2004), which, as we have checked, are 
consistent with those measured from best-fit continuum 
model spectra. The stellar mass estimates are total stellar 
masses derived fro m po pulation synthesis fits using the 
iBruzual & Charlotl (120031) models to SPSS bro a d-ban d 
photometry dKauffmann et alJ l2003bt ISalim et alJ 120071) . 
The adopted stellar mass estimate from photometry fits 
is a good indicator of the dynamical mass inside the ef- 
fective radius (Md yn oc o\R e lG, where R P is 
effective radiu s and G is 

Af* > 10'°M^ dBrinchmann & Ellisl l^OOl; iDrorv e taill2004t 
Padmanabhan et al. 2004); At the lower mass end of our 
sample (M» ~ 1O 9 5 M ), M* is ~ 0.2 dex larger than M d 
dDrory et al.ll2004h . 
To characterize how the effects of tidal interactions on host 



2 R e /G, where R e is the galaxy 
the gravitational constant) for 



yn 



http://www.mpa-garching.mpg.de/SDSS/DR7/ 



star formation and BH accretion evolve in AGN pairs as the 
merger progresses, we examine their statistical properties as a 
function of pair separation and velocity offset, and compare to 
control samples of single AGNs. We draw control AGNs from 
our parent sample of 138,070 AGNs. We match the control 
sample to the AGN pair sample in redshift to mitigate selec- 
tion biases in our flux-limited sample, and aperture bias due to 
the difference in the physical scales covered by the SDSS 3" 
fibers for galaxies at different redshifts. We also match stellar 
mass distribution to co ntrol mass-dependent ef fects of AGN 
host-galaxy properties (Kauffmann et al. 2003a). We show in 
Figure[T]the redshift and stellar mass distributions of the pair 
(tidal) sample and its control sample, respectively. We match 
the whole distribution rather than object by object. The con- 
trol samples are drawn to have the same distribution in red- 
shift and in stellar mass as the pair (tidal) sample by requiring 
that their Kolmogorov-Smirnov (KS) probabilities are larger 
than 50% at least. We have matched galaxy mass rather than 
luminosity in our control sample, because interacting galaxies 
often show younger stellar populations than do isolated galax- 
ies, maki ng galaxy lum inosity a biased indicator of mass (see 
also e.g jElhson et al.ll2008l) . 

To study the intrinsic dependence of host star formation and 
BH accretion on r p and Av in AGN pairs, we first verify in 
Figure |2]that host-galaxy stellar mass does not correlate with 
r p or Av, obviating the need to match control samples for each 
subset of AGN pairs at any fixed r p or Av. Figure|2]also shows 
that galaxy concentration, a measure of the mass distribution 
within galaxies, does not correlate with r p or Av, which is 
relevant because the distribution of mass within galaxies may 
also regulate merger timescales 6 . The r-band concentration 
index C r is defined as the ratio between the r-band Petrosian 

6 The mass distribution of dark matter halos at large radius dictates merger 
timescales on tens-of-kpc scales. 
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FIG. 3 . — Spectral indices H<5a and £>„(4000) as functions of interaction parameters for the pair (left block) and tidal (right block) samples. Symbols and line 
types are the same as in Figure [2] 



half-light radius and the radius e nclosing 90% of th e r-band 
Petrosian luminosity of a galaxy (Stra uss et al.l l2002). 

3. RESULTS: BLACK HOLE ACCRETION AND HOST-GALAXY 
RECENT STAR FORMATION IN AGN PAIRS 

In 33. II we present the distributions of BH accretion and 
host recent star formation properties of the pair and tidal sam- 
ples and their dependence on r p and Av, and compare them 
with control AGN samples. We show how these results de- 
pend on host-galaxy stellar mass ratio and concentration in 
jj3.2| § |3.3| focuses on the correlations between the interacting 
components in an AGN pair. We then examine in 33.4| the cor- 
relation between BH accretion and host recent star formation 
in AGN pairs and compare with those of control samples of 
single AGNs. 

3.1. Dependence on Projected Pair Separation and 
Recessional Velocity Offset 

3.1.1. Recent Star Formation in the Host Galaxies 

In this section, we address how galaxy tidal encounters af- 
fect recent star formation in AGN pairs. We cannot directly 
use the commonly used optical emission lines (Ha and [On]) 
as indicators of current star formation for the galaxies in our 
sample because they contain contributio n from AGN (e.g.. 
Ho et alJ [l997t iKauffmann etail l2003at iBrinchmann eTS] 



2004). Instead we adopt the continuum spectral indices H<5a 



and Z) K (4000) as indi c ators o f recent star formation, follow- 
ing IKauffmann et alJ d200 3a). Strong H5 absorption lines 
arise in galaxies that experienced a burst of star formation 
that ended ~0. 1—1 Gyr ago, and lar ge HSa values indi- 
cate strong recent starburst activ ity (e.g., Wor thev & Ottavianil 
[l997t IKauffmann et al]|2003bT) . D„(4000) measures the con- 
tinuum break around rest-frame 4000 A, which arises from a 
series of metal lines typical of 1-2-Gyr old stars (e .g. jBruzuall 
fl98llBalogh et all 1 9991: iBruzual & Charlotl2003h . £»„(4000) 
is small for young stellar populations and large for old, metal- 
rich galaxies, and can be u sed to characterize the l uminosity- 
weighted mean stellar age (Kauff mann et al.|[2003bl) . 

Figure[3]shows H<5a and D„(4000) as functions of r p and Av 
for the pair and tidal samples. Here and throughout, the error 
bars denote 1-er uncertainties in the medians. For AGN pairs 
with tidal features, both H<5 A and D„(4000) are significantly 



correlated with r p , despite having a large scatter. Spearman 
correlation tests show that the probability of null correlation 
P mU between H5 A (D„(4000)) and r p is 9 x 10" 5 (7 x 10" 6 ). 
The correlation probabilities quoted are estimated from the in- 
dividual data points directly. Pairs with smaller projected sep- 
arations on average exhibit stronger recent starburst activity 
and younger mean stellar ages. Compared to the control AGN 
sample matched in redshift and stellar mass distributions, 
AGN pairs with tidal features on average have larger H<5a and 
D„(4000) (median values of 1.4 ±0.1 and 1.55 ±0.01 com- 
pared to 0.1 ±0.1 and 1.67±0.01). The enhancement be- 
comes significant (> 3<r above the control) for objects with 
r p <10-30 h^ kpc (i.e., log(r p /h^ kpc) < 1.0-1.5). In 
a sample of ~ 10 s star-forming galaxies from SPS S DR4 
dAdelman-McCarthv et ail 12006b . IBrinchmann et al.1 (|2004|) 
find that D„(4000) anti-correlates with specific star forma- 
tion rate (SSFR) measured from optical emission lines. Ob- 
jects with the smallest separation in our sample (r p ~ 5 
kpc) have median D„(4000) — 1.40 ±0.04, correspond- 
ing to log (SFR/M H ,)/(yr~ 1 ) ~ -10.1 according to the cali- 
bration of IBrinchmann et al.1 (120041) . This is - 0.9 ± 0.2 
dex higher 7 than the SSFRs inferred for the control sam- 
ple (with median D„(4000) ~ 1.67 ±0.01, corresponding to 
log(SFR/M*)/(yr -1 ) of -11.0). 

We detect a similar dependence of H5a (D„(4000)) on r p 
in the parent population of all AGN pairs (i.e., regardless of 
the detection of tidal features). This suggests that the r p - 
dependence of H<5 A (D„(4000)) is robust against biases from 
the requirement of tidal features. However, the difference be- 
tween the median and overall H<5a (Ai(4000)) distributions 
of the pair and control samples is much smaller than that be- 
tween the tidal and control samples, most likely due to con- 
tamination from large-separation AGN pairs that are not yet 
interacting. On the other hand, our results seem to suggest that 
AGN pairs with tidal features show enhanced recent star for- 
mation at separations as large as ~ 100 h^' kpc. Another pos- 
sibility for this apparent enhancement at large separations is 
that AGN pairs with tidal features by construction are biased 
toward interacting systems with recent or more pronounced 

7 This enhancement is reduced to ~ 0.7 ± 0.2 dex if we adopt the cali- 
bratio n using mean instead of median values (Figure 1 1 of Brinchmann et al. 
2004). 
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FIG. 4. — Eddington ratio and [O III] A50O7 emission-line luminosity(in solar units) as functions of interaction parameters for the pair (left block) and tidal 
(right block) samples. Symbols and line types are the same as in Figure|2] 



star formation because younger stellar streams have higher 
surface brightness. 

In the tidal sample, we find no significant dependence of 
H(5a (Aj(4000)) on Av, after excluding small values of Av 
(< 30 km s" 1 ) which are dominated by redshift measurement 
uncertainties. Spearman correlation tests show that the prob- 
ability of null correlation between HSa (D„(4000)) and Av is 
0.72 (0.25). In the pair sample, we detect a weak dependence 
(Pnuii ~ 10" 3 ) of HS A (A,(4000)) on Av, such that pairs with 
smaller velocity offsets exhibit stronger recent starburst activ- 
ity and younger mean stellar ages. 

A caveat is that HS\ and D„(4000) are measured based on 
SDSS fiber spectra, which are subject to aperture bias due 
to the radial variation o f star formation within galaxies (e.g., 
iBrinchmann et al.ll2004l) . At the smallest and largest redshifts 
of our sample (z = 0.02 and 0.16), the 3" fiber corresponds 
to 1.2 kpc and 8.3 kpc, respectively, in the assumed cosmol- 
ogy. However, the r p distribution at any fixed redshift after 
correction for spectroscopic incompleteness does not depend 
significantly on redshift over the range considered, and there- 
fore our results on the r p dependence are robust against aper- 
ture bias. Repeating our analysis using SSFR es timates cor- 
rected for aperture bias (Brinchmann et al. 2004) provided in 
the MPA-JHU DR7 data product gives results consistent with 
those above 8 . 

We must caution that our results are limited to r p > 5 h^,, 
kpc. Our sample is incomplete at r p < 5 h^ kpc (see Paper 
I for details), because pairs with angular separations smaller 
than ~ 2" will not be resolved by the deblending algorithm 
of SDSS photometry and therefore will not be included in our 
sample; at z ^0.16, i.e., the highest redshift bin of our sample, 
this corresponds to ~ 5 h^' kpc. We estimate that the incom- 
pleteness level due to this resolution limit at r p = 5 h^, kpc 
is ~ 1-5% (due to the highest redshift bin around z = 0.16), 
according to the number of smaller separations pairs (r p ~ 1- 
5 kpc) resolved at lower redshift in our sample assuming 

8 The correction is made by applying the SSFR distrib ution at any given 
j> — r and r-i color to SDSS photometry outside the fiber (Brinchmann et al. 
12001 . We use these estimates for verification purposes only, and take the di- 
rectly measured spectral indices (after subtracting emission lines) as our base- 
line values becau se they are less prone to population-dependent bias (e.g., 
ISalim et aU200l) . 



that the r p distribution around r„ 



h 7 Q kpc does not vary 



with redshift. This incompleteness level increases with de- 
creasing r p , and it reaches ~ 100% at r p ~ 0.8 hyg kpc (i.e.. 



the physical scale corresponding to < 
redshift boundary of our sample). 



■ 2" at z = 0.02, the lower 



3.1.2. Black Hole Accretion 

We now address how galaxy tidal interactions affect BH ac- 
cretion in the early stage of mergers in AGN pairs, and com- 
pare these effects with that of recent star formation to con- 
strain their causal relationship. We adopt the [O in] A5007 
emission-line luminosity L[om] and the inferred Eddington ra- 
tio (LBoi/^Edd) as indicators of BH accretion activity. We cal- 
culate L[oni] based on fluxes provided in the MPA-JHU cata- 
logue from Gaussian fits of continuum-subtracted emission- 
line spectra. The forbidden line [O in] A5007 arises from 
the narrow-line region 9 , and its luminosity is correlated with 
broad-ban d continuum luminosities in unobscured AGNs 
(e.g.. iKauffmann et alj l2003at iReves etal] 120081: IShen et alj 
2011b), albeit with a si gnificant scatter (e. g., 0.36 dex in the 
correla tion with M2500; Rey es et all l2008h . iHeckman et alJ 
d2005l) found that [Om]A5007 luminosity correlates well 
with the 3-20 keV luminosity with a dynamic range in lu- 
minosity of 4 dex in a sample of 47 hard X-ray selected lo- 
cal AGNs. To infer bolometric AGN luminosity LboI, we 
use L[oni] uncorrected for dust extincti on and the appropri- 
ate bolometric correction factor (3 500; IHeckman et al.ll2005l : 
iLiu et al.ll2009l: IShen et alJl201 lbl) . We do not correct L [0 ni] 
for dust extinction from the hosts in view of the uncertainties 
in the assumed reddening laws, which rely o n over-simplified 
dust-screen models (e.g., Reyes et al. 2008). For the Edding- 
ton luminosity, we adopt a BH mass (Mbh) inferred from the 



We have subtracted the contribution from star formation to Li 



•I o 11 1 ] 

based on the position o n the BPT diagram (Baldwin et a i] |l981t 
Veilleux & Osterbrock 1987) which features the diagnostic emission-line ra- 
tios [O III] A5007/H/3 and [N II] A6584/Ha, following the empirical method 
of Kauffmann & Heckman 1 2009). This approach constructs composite ob- 
jects using a linear combination of the spectra of pure star-forming galaxies 
and pure AGNs. Templates of pure star-forming galaxies are made with ob- 
jects located at the lower edge of the star-forming sequence, whereas pure 
AGNs are built using those located at the farthest ends of the Seyfert and 
LINER regimes. Trajectories on the diagnostic diagra m are th en produced 
from the constructed composites (Kauffmann & Heckman 2009). 
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stellar velocity dispersion (er*) measured from the SDSS spec- 
trum, assuming th e Mbh-c* relation observed in local inac- 
tive g alaxies (e.g., iFerrarese & MerrittllgOOOt iGebhardt et all 
120001) and calibrated bv lTremaine et al.1 (120021) . As we will 
discuss below, this is likely an overestimate of the true mass 
for small-separation AGN pairs due to the effects of tidal in- 
teractions on the internal kinematics of the host stellar popu- 
lations. 

Figure |4] shows [Oiii] A5007 luminosity and Eddington ra- 
tio as functions of r p and Av. In the tidal sample, the strength 
of BH accretion activity increases with decreasing r p . There is 
also a similar dependence seen with Av after excluding small 
values (< 30 km s" 1 ) which are dominated by redshift mea- 
surement uncertainties. Spearman correlation tests show that 
the probability of null correlation between L[oni] (iBoi/^Edd) 
and r p is 10~ 3 (2 x 10~ 3 ). These weak correlations are not 
as significant as those found between recent star formation 
indicators and r p . The enhancement in AGN activity is sig- 
nificant for objects with r p < 10-30 h™ kpc (i.e., log^/h^ 
kpc) < 1.0-1.5). Objects with the smallest separation in 
our sample (r p ~ 5 hy, 1 , kpc) have median log(L[om]/Zo) and 
log(L Bo i/iEdd) values ~ 0.7 ±0.1 and ~ 0.5 ±0.1 dex larger 
than the median values of the control AGNs. AGN pairs with 
tidal features on average exhibit accretion activity stronger by 
~0.2-0.3 dex than control AGNs matched in redshift and stel- 
lar mass distributions. 

In the pair sample we detect a similar dependence of L[om] 
(LBoi/^Edd) on r p- This verifies that the r p -dependence re- 
sult is robust against biases from our identification of tidal 
features. However, the differences in the median value and 
distribution of Lpm] (^Boi/^Edd) between the pair and control 
samples are much smaller than those between the tidal and 
control samples. The dependences of L[om] and LBoi/^Edd on 
r p are less prominent than that of recent star formation shown 
in Figure [3] and therefore can be more difficult to detect. 

Figure [5] shows that er„ at small r p (~ 10-30 h^' kpc) is 
on average larger than that of large separation pairs and the 
median value of control AGNs by ~ 20 ± 5 km s" 1 in both 
the pair and tidal samples. The stellar mass and effective ra- 
dius are independent of r p . Therefore the larger cr* cannot 
be due to deeper potentials; it is most likely the effect of the 
interactions with the transformation of orbital into internal en- 
ergy. As a result, the inferred BH mass is likely overestimated 
by ~ 0.2 dex on average, and by extension, the enhancement 
of iBoi/^Edd with decreasing r p is likely underestimated by 
~ 0.2 dex in AGN pairs with the smallest separations in our 
samples. We see a similar enhancement in gas velocity disper- 
sion (measured from Gaussian fits to the forbidden emission 
lines) as well (Figure |5). 

Because the redshift distribution corrected for spectro- 
scopic incompleteness at any fixed r p does not depend signif- 
icantly on r p (Paper I), our results involving luminosities are 
robust against redshift-induced biases. Repeating our analysis 
using [Oiii] A5007 equivalent width (EW) as an indicator of 
BH accretion yields consistent results. Another caveat is that 
the [Oiii] A5007 fluxes are based on SDSS fiber spectra and 
could therefore be subject to aperture bias (see the discussion 
at the end of $3.1. 11 1. The typical size of the narrow-line re- 
gions of lower-luminosity AGNs is several hundred parsecs, 
which is smaller than the aperture size of objects with the low- 
est redshift in the sample, and therefore the bulk of narrow- 
line region [O in] A5007 emission should be included in the 
fiber spectra of all objects in the sample. In addition, we have 



verified that the distribution of separations for 5 < r p < 100 
hyy kpc at each redshift bin (after correction for spectro- 
scopic incompleteness) does not vary significantly with red- 
shift (over the range that we consider, i.e., 0.02 < z < 0.16). 
Therefore, the correlations of the galaxy properties with r p are 
not affected by redshift dependent biases. 

3.2. Dependence on Host-Galaxy Mass and Structure 

We now address how host-galaxy internal properties regu- 
late the effects of tidal interactions on recent star formation 
and AGN activity in AGN pairs. We examine the role of mass 
ratio in §13.2.11 and that of gal axy structure (as indicated by 
concentration index) in $3.2.21 We focus on the dependence 
of recent star formation and BH accretion as a function of r p . 

3.2.1. Effect of Host-Galaxy Mass Ratio: Major versus Minor 

Mergers 

Mergers and close interactions (hereafter mergers for short) 
of galaxy pairs with comparable masses (so called "major 
mergers") as a driver for star formation and AGN activ- 
ity have been the topic of extensive th eoretical and obser- 
vational research (e .g., see reviews by iBarnes & Her nquist 
Il992t IStruckl [2006). It is generally accepted that major 
mergers of gas-rich spirals are responsible for the most 
luminous starburst galaxies , including the ne a rby ultra- 
luminous infrared galaxies (ISanders et al.l 119881 iHernquistl 
119891: iMihos & Hernquistll 19961) and s ub-millimeter galaxies 
at higher redshifts (Blain et alT l2002b IConselice etaLl l2003t 
iTacconi et al.ll2008l) . Minor mergers are more frequent than 
major mergers. At z < 1 there is evidence that minor mergers 
may play a dominant role i n inducing moderat e star forma- 
tion in early-type galaxies (Kavir ajet al.ll201lh and in driv- 
ing the structural evoluti on of intermediate-mass galaxies 
(Lopez-Saniuan etai1l2010l) . It is still debated whether ma- 
jor or minor mergers or internal instabilities are the p rimary 
trigge rs for various AGN populations (for a review see Jogee 
2006). 

We divide the our AGN pair samples into major and minor 
mergers, i.e., with stellar mass ratios smaller and larger than 
3:1, respectively. The stellar mass ratios of our sample span 
the range 1 < M*,i /M* 2 < 10, where M„.i denotes the stellar 
mass of the primary and M* 2 is that of the secondary galaxy. 
38% of our AGN pairs have mass ratios smaller than 3. The 
sample is incomplete for pairs with mass ratios larger than 
~10, because of the limited dynamical range of the parent 
spectroscopic galaxy sample. 

Figure [6] shows the recent star formation and AGN indica- 
tors as functions of r p for major and minor mergers for the pair 
and tidal samples. In the pair sample, galaxies in major merg- 
ers have median stellar masses of ~ 10 I0 ' 7 M Q , whereas the 
primary and secondary of minor mergers have median stel- 
lar masses of - 10 10 ' 2 M Q and - 10 110 M Q , respectively. We 
observe a dependence of host recent star formation and AGN 
activity on r p in both major and minor mergers. In minor 
mergers, we detect and enhancement of recent star formation 
and AGN activity indicators at small r p in both the primary 
and secondary galaxies. AGN pairs with separations of ~ 5 
h^Q kpc have median SSFRs that are larger than those with 
separations of r p ~ 100 h^,' kpc (approaching the case of con- 
trol samples of single AGNs) by ~ 0.8, ~ 0.6, and ~ 0.7 
dex (with typical uncertainties of 0.2 dex) in major, minor- 
primary, and minor-secondary gala xies respectively, based o n 
the SSFR-D„(4000) calibration of BrinchmamTetalJ |2004). 
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FIG. 5. — Stellar and gas velocity dispersions as functions of interaction parameters for the pair (left block) and tidal (right block) samples. Symbols and line 
types are the same as in Figure [2] 



The enhancement level of median L[oni] (^-Boi/^Edd) m AGN 
pairs from r p ~ 100 to ~ 5 h^, 1 , kpc is 0.4, 0.5, and 1.2 dex 
(0.4, 0.6, and 0.7 dex) for major, minor-primary, and minor- 
secondary galaxies, respectively, with typical uncertainties of 
0.2, 0.4, and 0.3 dex. The most notable difference between the 
secondary and primary in minor mergers is that the enhance- 
ment level of median L[oni] from r p ~ 100 to 



5 kpc 



is significantly higher in the secondary than in the primary 
galaxy. The difference between the secondary and primary in 
the enhancement of LBoi/^Edd is less significant, although we 
find that the enhancement of cr* is larger in the secondary than 
in the primary which may cause LBoi/^Edd to be more severely 
underestimated (see the discussion at the end of ^3. 1.2b . In the 
tidal sample, we find similar results to those in the pair sam- 
ple, except that there is a significantly larger enhancement of 
D„(4000) in the secondary than in the primary of minor merg- 
ers. We discuss the implications of these results in §H] 

3.2.2. Effect of Host-Galaxy Structure: Bulge versus Disk Mergers 

The distribution of mass within galaxies regulates when 
gravitational instabilities are induced in galaxy tidal interac- 
tions. Simulations of major mergers of disk galaxies show 
that a dynamically hot bulge component stabilizes a disk 
against bar instability, thu s delaying merger-induced activity 
dMihos & Hernquistlll996l) . We address whether such an ef- 
fect is present in AGN pairs. We use the r-band concentration 
index, C r , as a measure of mass distribution. IStrateva et al.l 
d2001l) have shown that C r correlates with ga laxy morpholog- 
ical ty pe albeit with a large scatter (see also IShimasaku et all 
2001); elliptical galaxies have larger C,- than do spirals. C ~ 
5.5 for a pure de Vaucouleurs elliptical whereas C ~ 2.3 for 
a pure exponential disk. In particular, Stra teva et al.l (1200 ll) 
have found that a value of C r = 2.83 divides early- and late- 
type galaxies with a reliability of 83% for each and with com- 
pleteness of 70% and 80%, respectively. As shown in Figure 
[7] the C r distribution of our AGN pair samples is close to be- 
ing a Gaussian with a median value of C r = 2.8, and lies in 
the intermediate region between early- and late-type galaxies. 
We divide the AGN pair galaxies into those with C r > 2.8 and 
those with C r < 2.8 (each component in a pair is treated sep- 
arately). 

We first examine how host-galaxy structure regulates the ef- 



fects of tidal interactions on host star formation in AGN pairs. 
Figure [7] shows the dependence of recent star formation and 
AGN activity indicators on r p for early- and late-type galaxies 
in the pair and tidal samples, respectively. In both samples we 
find a clear dependence of host star formation indicators on 
r p in AGNs with late-type host galaxies. From r p ~ 100 to 
~ 5 hjQ kpc, the median SSFR values are enhanced by ~ 0.4 
dex in AGNs with late-t ype hosts, according to the SSFR- 
D„(4000) calibration of iBrinchmann etal.1 (|2004| ). In AGNs 
with early-type hosts, however, the trend is much less promi- 
nent. In the pair sample, neither H5\ and D„(4000) show 
much dependence on r p until r p < 10 h^ kpc, whereas the 
dependence becomes prominent starting at r p < 30 hyg kpc 
in late-type hosts. In the tidal sample, there is little depen- 
dence of H(5a and D„(4000) on r p . Our comparison between 
the early- and late-type hosts of AGN pairs suggest that tidally 
enhanced star formation happens earlier and more prominent 
in late- than in early-type hosts of AGN pairs. 

We now assess how host-galaxy structure regulates the ef- 
fects of tidal interactions on BH accretion in AGN pairs. As 
shown in Figure [7] the median L[oni] (Lboi/Leu) value in- 
creases with decreasing r p in AGNs with both early- and 
late-type host galaxies in both the pair and tidal samples. 
In the pair sample, the median Lpm] value increases from 
r p - 100 to - 5 h"^ kpc by 0.7 ±0.1 dex (1.0±0.1 dex) of 
AGNs with early-type (late-type) hosts, respectively (here we 
neglect the smallest r p bin in early-type hosts to focus on the 
overall trend). The increase of the median LBoi/^Edd value 
seems less prominent: 0.7±0.1 dex (0.4±0.1 dex) for AGNs 
in early-type (late-type) hosts. However, the enhancement of 
i-Boi/i-Edd could be underestimated because cr*, and by exten- 
sion Mbh, is likely overestimated for AGN pairs with small 
separations (Figure|5]i; late-type hosts show more <r* enhance- 
ment. We detect similar trends in the tidal sample. The differ- 
ent levels of L[oni] median enhancement in AGNs with early- 
and late-type hosts suggest that AGNs with late-type hosts are 
more subject to tidally enhanced BH accretion than those with 
early-type hosts. 

We have verified that the trend we see in the early and late- 
type hosts of AGN pairs in indeed due to the difference in their 
concentration rather than a mass selection effect. First, we 
compared the stellar masses of the C r > 2.8 and the C r < 2.8 



8 



Liu, Shen, & Strauss 




tidal" i-o 

0.5 



0.0 




i~L , , H o.o E 




10 11 12 10 11 12 

major" iog(M./M ) "minor" log(M,/M e ) 



10 11 12 10 11 12 

major" log(M./M s ) "minor" i°g(M./M Q ) 



2.0 



o 
o 

° 1.6 



1.2 



1. 5 



^ -3.0 
o 




cn 
O 



















+H+ 


-H-H- 


+ftf; 






-T-H+ 


-++++- 
























-++++- 


-r-r-r-r- 




+ 




-r-r-r-r- 


-++++- 




2.0 



o 
o 

° 1.6 



1.2 



1 .5 



^ -3.0 

o 



CP 

o 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 







+H+ 










-++++- 


-++++- 




-fc 




1 1 1 1 [ 1 1 1 1 


4+++- 


















+T++- 


+T-T-T- 


+T-T-T- 




+ 


d+T-T- 


-++++- 



_L 



I I I I I I I I I 



















+T++- 


+T-T-T- 


+T-T-T- 


-F 




-++++- 


-++++- 


++++- 



_L 



i i i i I i i i i 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 
'ogOp/hyo k Pc) log(r p /h^ kpc) 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 
'ogOp/h^ k Pc) log(r p /h^ kpc) 



FIG. 6. — The effect of host-galaxy mass ratio on tidally enhanced star formation and AGN activity. The left (right) block is for the pair (tidal) sample. In each 
block, the left column is for major mergers with mass ratios smaller than 3, whereas the right column is for minor mergers with mass ratios larger than 3. For 
minor mergers, we show the primary (open circles in magenta) and secondary (open triangles in green) galaxies separately. Small dots denote individual objects, 
whereas large symbols represent median values at each r p bin. Error bars indicate uncertainties of the median values. The panels in the first row display stellar 
mass distributions. 



samples (here we quote the numbers for the pair sample) and 
found that they defer by only 0.2 dex on average; their stellar 
mass distributions are consistent with being Gaussian distri- 
butions, with the C r > 2.8 sample having a median mass of 
10.81 (in units of log(M*/M©)) and a standard deviation of 
0.34, and the C r < 2.8 sample having a median stellar mass 
of 10.62 and a standard deviation of 0.43. This 0.2 dex differ- 
ence is much smaller than that between the primaries and sec- 
ondaries (~ 1 dex) in minor merger pairs (§ 13.2.1) . While the 
C r > 2.8 sample has on average larger stellar masses than the 
C r < 2.8 sample (as expected due to the mass-concentration 
correlation), the stellar mass difference seems to be too small 
to be the dominant factor of the trend we see. Second, to test 
whether the effect of this small mass difference is indeed mi- 
nor, we redid our analysis using subsets of the C r > 2.8 and 
the C r < 2.8 samples matched in stellar mass distribution. The 
trend was still present, verifying that it was indeed due to con- 



centration effect rather than a mass selection effect. 

3.3. Correlations between the Interacting Components in an 

AGN Pair 

In this section we examine the correlations between the two 
interacting components in each AGN pair. The statistical sam- 
ple of AGN pairs presented in Paper I allows us to address 
whether there is a correlation in the strength of BH accre- 
tion between the interacting components. We show in Figure 
[8] the correlations in indicators of recent star formation be- 
tween the interacting components for the pair and tidal sam- 
ples. We plot H<5a (D„(4000)) of one component (randomly 
pulled from the two) against that of the other component in an 
interacting pair. Also shown for comparison are density con- 
tours from control AGNs. We draw control samples of AGNs 
for both components, each of which is matched in redshift and 
stellar mass distributions to one component of the AGN pairs. 
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FIG. 7. — The effect of host-galaxy structure on tidally enhanced star formation and AGN activity. The left (right) block is for the pair (tidal) sample. In each 
block, the left column is for galaxies with r-band concentration index C r > 2.8 (corresponding roughly to ellipticals) whereas the right column is for those with 
C r < 2.8 (spirals). Small dots denote individual objects, whereas large symbols represent median values at each r p bin. Error bars indicate uncertainties of the 
median values. The panels in the first row show the C,- distributions. 



The two control samples are then randomly "paired" with the 
requirement that each control "pair" has a LOS velocity offset 
smaller than 600 km s" 1 , although our results remain the same 
without this constraint. As shown in Figure [8] there is no cor- 
relation between the two components in the control sample; 
Spearman null probabilities are close to 1. We have run 1000 
Monte Carlo simulations for the control sample and verified 
that the Spearman correlation probabilities shown in Figure 
[8] are typical. For the pair and tidal samples, on the other 
hand, we detect a statistically significant correlation between 
the recent star formation indicators of the interacting galax- 
ies. Spearman tests show the probability of null correlation 
is 1 x 10" 4 (1 x 10" 4 ) in H5 A (£>„(4000)) for the pair sam- 
ple; the null probability is even smaller for the tidal sample, 
2 x 10" 5 (5 x 10" 5 ) in HS A (D„(4000)). There is a substantial 
scatter in the correlations though. The correlations seen in 
Figure [8] for AGN pairs are reminiscent of the color correla- 



tion s een between inactive galaxy pai r s (the Holmberg effec t; 
e.g., lHolmberg|[T958l lMadordfl986l iKennicutt 



aberg errect; 
etafl[l987l) . 



which is generally a ccepted as evidence for me rger-induced 
star formation (e.g., Laurikain en & Moleslll989l) . 

We now address whether there is an effect in BH accretion 
activity similar to that observed in recent star formation. To 
avoid selection bias due to the common distance moduli in 
a pair in a flux-limited sample, here we adopt the rest-frame 
[O in] A5007 equivalent width (EW) as an indicator of BH ac- 
cretion, instead of [O in] A5007 luminosity or Eddington ra- 
tio. Figure [9] shows the [O in] A5007 EW of one component 
(randomly pulled from each AGN pair) as a function of that 
of the other component in the pair for the pair and tidal sam- 
ples. In the pair sample, we detect a statistically significant 
correlation in the [O in] A5007 EW, although the scatter is sig- 
nificant (~ 1 dex); the Spearman null probability is 4 x 10~ 5 . 
The correlation in the tidal sample is even stronger, with a 
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FIG. 9. — As in Figure [8] but for the correlation of the rest-frame 
[O III] A5007 EW (an indicator of BH accretion activity) between the in- 
teracting components in each AGN pair. The left (right) panel is for the pair 
(tidal) sample. 

null probability of 5 x 10~ 6 . We find no correlation between 
[O in] A5007 EWs of the two components in the control sam- 
ples; the Spearman null probabilities are close to 1. The cor- 
relation of AGN activity between the two components in an 
AGN pair suggests that there is tidally enhanced AGN activ- 
ity in both components. Alternatively, a single active SMBH 
ionizing both galaxies could also produce such a correlation, 
at least for the subset of AGN pairs with separations smaller 
than galaxy sizes. The latter scenario is unlikely, however, as 
we find no significant correlation 10 in the diagnostic emission- 
line ratios [Om]/H/3 11 between the interacting components 
as shown in Figure [To] suggesting that the two components 

10 There is a tentative correlation of [N II]/Ha in both the pair and tidal 
samples, which could be due to the tidal ef fect on the ga s-phase metallicity 
observed in inactive galaxy pairs (e.g., Ellison et al. 2008). 

1 1 The [O [O IIIJ/H/3 ratios have not been corrected for contributions from 
star formation. The results show no significant change when [O III] is cor- 
rected for star formation. 



are not powered by the same central engine. 

3.4. Correlations between Recent Star Formation and Black 
Hole Accretion 

AGN luminosity is correlated with H5a and D„(4000) in 
the host galaxies of z ~ 0.1 AGNs, in the sense that the 
hosts of higher-luminosity AGNs have younger stellar pop- 
ulatio ns and higher starburst fractions in the last 0.1 -1 Gyr 
(e.g.. lKauffmann et alJl2003at lNetzerll2009l: lLidl2010l) . Here 
we address whether AGN luminosity is correlated with recent 
star formation in hosts in AGN pairs, and how such a correla- 
tion, if present, compares to that in ordinary AGNs. 

FigurefTTIdisplays L[omj andLBoi/^Edd as a function of H<5a 
and D„(4000) for the pair and tidal samples. Also shown for 
comparison are contours from control AGN samples matched 
in redshift and stellar mass distributions. There is a cor- 
relation between Lroni] and H<5a (A,(4000)) when H(5a> 2 
(D„(4000)< 1.6) for the pair sample and when US A > -2 
(D„(4000)< 1.9) for the tidal sample, respectively. Lfioi/^Edd 
is also correlated with H<5a and D„(4000) for both samples, 
although at least part of the correlation is driven by the mass 
dependence of both er* (and by extension Lndd) and H<5a 
(D„(4000)). The relation between AGN luminosity and H<5 A 
(D„(4000)) of the pair sample is almost identical to that of 
the control AGN sample. While the tidal sample occupies 
a similar scaling relation between AGN luminosity and H<5a 
(D„(4000)) to that of the control sample, it is skewed to- 
wards higher AGN luminosities and larger starburst fractions 
(younger mean stellar ages). 

4. DISCUSSION 

We discuss implications of our results on tidally enhanced 
star formation and BH accretion. We compare our results with 
previous observations of inactive galaxy pairs and of single 
AGNs in galaxy pairs in $4.11 and to predictions from galaxy 
merger simulations in $4.21 

4.1. Comparison with Observations of Galaxy Interactions 

4.1.1. Tidally Enhanced Star Formation 

We have detected a correlation between H8a (and 
D„(4000)) and projected separation r p in AGN pairs; sys- 
tems with smaller separations have larger H<5a and smaller 
D„(4000) indicative of higher SSFRs, with the enhancement 
becoming significant for r p < 10-30 h^' kpc (Figure [3}. We 
also find a weak correlation between H<5a (and D„(4000)) 
and Av in both our tidal and pair samples, after exclud- 
ing small values of Av which are dominated by redshift 
measurement errors. These results are in broad agreement 
with previous findi ngs based on statistical samples of in- 
active galaxy pairs ([Barton et al. 2000; Lambas et al.l 120031: 
Alonso et all 120041; INikolic et all 12001; lEllison etafl 120081; 
Li et al . 2008a; Darg et a02010al) . For example, in a volume- 
limited sample of 12,4 92 main galaxies from the SDSS DR1 
( Abazaiian et alj 120031) selected to have ph otometric com- 
panions within 300 kpc, INikolic et all d2004l) used SFRs in- 
ferred from Ha luminosities to find that the mean SSFR is 
significantly enhanced for r p < 30 h^ kpc. These authors 
also find a weak anti-correlation betw een SSFR and reces - 
sional velo city difference Av (see also Lambas et al. 2003). 
Similarly, Ellis on et al.l (|2008) found an enhancement in the 
star formation rates of galaxy pairs at r p < 30^40 h^,' kpc 
in a study of 1716 emission-line galaxies se lected from the 
SDSS DR4 (Adelman-M cCarthv et ai]|2006l) with compan- 
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FIG. 10. — As in Figure|8] but for the correlation of the diagnostic emission- 
line ratios between the interacting components in each AGN pair. The top 
(bottom) row is for the pair (tidal) sample. 

ions within Av < 500 km s" 1 , r p < 80 kpc, and stellar 
mass ratios smaller than 10. 

The average SSFR enhancement level we observe in AGN 
pairs is 0.7-0.9 ± 0.2 dex from r p > 100 to ~ 5 h^ kpc. This 
is comparable to some previous estimates based on inactive 
gala xy pairs (e.g., f actors of 1.5-4 from r p > 100 to r p = 20 
kpc. lLi etaT 2008a), but seems to be la rge r than others taken 
at face value (e.g.. lAlonso et af]|2004l and iDarg etail 1201 Oal 
found an enhancement of a factor of 2). The variation in the 
observed enhancement levels may be due to difference s in the 
galaxy samples studied. For example, the lAlonso et al.l(l2004l) 
results are based on galaxy pairs selected in galaxy groups or 
clusters with virial masses in the rage of 10 13 -10 15 M Q and 
therefore may be biased towards early-type galaxies. Indeed, 
we find tentative evidence that the star formation enhance- 
ment is more prominent and occurs at larger separations in 
late-type than in early-type galaxies (Figure [7}. The variation 
in the observed enhancement levels may also be due to dif- 
ferences in the interaction stages examined. For ex ample, the 
factor of 2 difference inferred by IDarg et aD (1201 Oal) refers to 
the average difference between their merger sample (includ- 
ing various separations hence interaction stages) and the con- 
trol sample, which may be much smaller than the difference 
observed for small separation pairs. Other factors that may 
affect the observed level of enhancement include how control 
samples are defined and the contamination from closely sepa- 
rated galaxy pairs that are not tidally interacting. 

We find that the enhancement level of recent star forma- 
tion in AGN pairs with host stellar mass ratios M^i/M*^ > 
3 (minor mergers or interactions) is comparable to (if not 
slightly smaller than) that in those with stellar mass ratios 
i /M* 2 < 3 (major mergers or interactions), at least when 
the two galaxies are separated by more than r p ~ 5 h^,' kpc. In 
addition, the relative enhancement level is similar in both the 
primaries and the secondaries in AGN pairs with ste llar mass 
ratios M*,. /M* j2 > 3. Similarly. [Ellison et al.l (120081) detected 
SFR enhancement in their sample of inactive galaxy pairs 
drawn from the SDSS whose stellar masses are within a factor 
of 10 of each other. On the other hand, based o n 167 galaxy 
pairs drawn from the CfA2 Redshift Survey dHuchra et al.l 



1990, 1995), Wood s et all d2006l) suggest that galaxies in ma- 
jor interactions are more likely to show enhanced star for- 
mation activity than those in minor interactions with abso- 
lute magnitude differences larger than 2. In a sample of 
1204 galaxies in spectroscopi c pairs and compact groups se - 
lected from the SDSS DR5 dAdelman-McCarthy et afll2007l) . 
IWoods & Gellerl d2007T) find that in minor mergers the sec- 
ondary, but not the primary shows tidally triggered star for- 
mation. On the other hand, in a study of 1258 spectroscopi c 
galaxy pairs d r awn fr om the 2dF survey (IColless et al.l! 200 1 ). 
Lamb aset al.l ((2003) find that in minor mergers the bright 
components are more likely to show tidally-enhanced star for- 
mation than the faint components. We speculate that these 
contradictory results could be due to difference in the galaxy 
samples studied, variation in their sample sizes, and criteria 
to define major and minor mergers. For exampl e, we adopt 
stella r mass rati o in the curren t study ( see al so Elliso net alj 
120081). whereas iLambas et al.1 d2003l) . IWoods et all d2006l) . 
and I Woods & Gellerl (12007b employed the ratio of galaxy lu- 
minosities in their studies. 

We caution that our conclusion, that recent star forma- 
tion is enhanced in AGN pairs during the early merger 
stages, may not hold for all types of galaxies, given the 
inherent limitations of our sample selection. Low-redshift 
optical AGNs reside in galaxies with stellar mass, galaxy 
structure, and star formation properties intermediate be- 
tween you ng disk-dominated and o ld bulge-dominated galax - 
ies (e.g., iKauffmann eTaT] l2003allcl : iHeckman etaD 120041) . 
Whether tidal encounters enhance star formation in galaxies 
with all ranges of galaxy mass and structure is beyond the 
scope of this paper, although there is mounting evidence that 
the specific levels of star formation enhance ment depend on 
interaction and/or progenitor properties (e .g., Bushouse ]! 19861 : 
Bushouse et al.ll988l:lBergvall et al.l2p03tlBrosch et alJ2004t 
Woods et al.ll2006l:lRogers et al.ll2009l) . 



4.1.2. Tidally Enhanced Black Hole Accretion 

Previous statistical studies of galaxy pairs found ev- 
idence for a higher fraction of AGNs in paired than 
in is olated galaxies dKennicutt et al J 119871: lAlonso et alJ 
1 20071: IWoods & Gellerl2007l:lRogers et al.ll2009blEilison et al.l 
1201 1|) where a s oth e rs detected no sign i ficant difference 
dBarton et al.1 120001: lEllison et all 120081: iLi et all l2008bt 
Pare et al.l 1201 Oal) . We have found a correlation between 
[O in] A5007 luminosity (and Eddington ratio) and projected 
separation in AGN pairs; systems with smaller separations 
host more powerful AGNs and the enhancement becomes sig- 
nificant for r p < 10-30 h^,' kpc (Figure |4j, although the en- 
hancement is only moderate, which is ~ 0.2-0.3 dex on aver- 
age in luminosity/Eddington ratio relative to control samples 
and rises to ~ 0.5-0.7 dex for the closest pairs in our sample 



5 h^ 1 , kpc). 



We find that the relative enhancement in 



AGN luminosity from r p > 100 to ~ 5 h 7 g kpc is higher in 
the secondary than in the primary in AGN pairs with stellar 
mass ratios M*,i/M».2 > 3. Such a difference may be driven 
by the fact that less massive galaxies on average have higher 
gas fractions, or that they are more prone to gravitational in- 
stabilities. The offset Seyfert 2 nu cleus in the mino r merger 
system NGC 3341 discovered by iBarth et all (T2008) may be 
a manifestation of the same effect, where the BH in the sec- 
ondary is experiencing much more accretion than that in the 
primary. 

We discuss several reasons why the causal link between 
galaxy interactions and AGN activity has been controversial. 
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indicate control samples of single AGNs matched in redshift and in stellar mass. 



First, the merger- AGN connection see ms to depend sensi- 
tively on AGN luminosity. For example, ISerber et al.l (120061) 
analyzed ~ 10 7 close neighbors of ~ 2,000 SDSS quasars at 
Z < 0.4 with median LboI ~ 10 45 5 erg s" 1 and found that they 
reside in higher local (< 100 kpc) over-density regions than 
do normal L* galaxies. In a much lower AGN lum inosity 
regime (median L Bol ~ 10 430 erg s'^ JLi et alJd2006b studied 
cross correlations of ~ 90,000 narr ow-line AGNs with galax- 
ies at z ~ 0.1 from the SDSS DR4 ( I Adelman-McCarthv etall 
2006), and found that on scales < 70 kpc, AGNs cluster 
marginally more strongly than control samples of inactive 
galaxies matched in redshift, stellar mass, concentration, ve- 
locity dispersion, and mean stellar age, whereas they cluster 
more weakly than inactive galaxies on 0.1-1 Mpc scales, and 
cluster similarly on scales larger than a few Mpc. The AGN 
luminos i ty reg ime we are probing here is similar to that of 
iLi et all (12006). As we show in Figure [12] the fraction of 
galaxy pairs hosting either single or double AGNs 12 among 
all galaxy pairs /agn shows no significant dependence on r p 
from ~ 100 h^ kpc to > 5 h^ kpc fo r all AGNs. T his is 
consistent with the clustering results of ILi et all (120061) . On 
the other hand, if we consider only the stronger AGNs (e.g., 
L[om] > 10 40 erg s" 1 corresponding to LboI > 10 43 5 erg s" 1 ), 
/agn rises prominently with decreasing r p for r p < 30 hjQ kpc, 
with a factor of ~ 2.5 ±0.5 (~ 10 ±2) increase from ~ 100 



kpc to > 5 h 7 y kpc for those hosting either single or dou- 



h7u 

ble (those hosting double) AGNs. This comparison suggests 
that the overall increase we have detected in AGN luminos- 
ity/Eddington ratio in close AGN pairs is a result of tidally 
enhanced BH accretion events on top of a background of weak 
AGN events driven by secular processes not associated with 
interactions. 

Second, the r p correlation we detect for AGN luminosity is 
not as prominent and significant as for indicators for recent 
star formation. Therefore we speculate that small samples of 
AGN pairs with tens of objects or even fewer may simply be 



12 To make a fair comparison with the Li et al. (2006) clustering results, 
we adopt the same AGN diagnostic criteria here, i.e., AGNs according to 
fhe lKauf fmann et al. i 2003a) empirical curve with all four emission lines de- 
tected with S/N> 3. 



insufficient to reveal such a weak correlation. In addition, the 
difference between the overall distributions of AGN indica- 
tors in the pair and control samples is smaller than that in the 
tidal sample (Figure|4]i. Thus, the effects of tidal interactions 
on AGN activity may be underestimated due to contamina- 
tion by closely separated galaxy pairs that are not undergoing 
tidal interactions. Furthermore, we have drawn control sam- 
ples with identical redshift and stellar mass distributions. Be- 
cause galaxy interactions also enhance star formation activity, 
their effects on AGN activity may be underestimated if con- 
trol samples are matched in c olor or in stella r age, as was done 
in the clustering analysis of ILi et all d2006). Finally, the en- 
hancement in AGN activity we observe becomes prominent 
only when studied as a function of separation and when r p 
is sufficiently small (< 10-30 hyQ kpc). As a result, the ef- 
fect will not be detectable in samples of galaxy pairs domi- 
nated b y larger-separation systems. For example, iDarg et all 
( 20 1 0a) find tentative evidence for a slightly higher AGN frac- 
tion in mergers identified for tidal features using SDSS images 
at later stages (32% at r p ~ 5 kpc) than at earlier stages (23% 
at r p ~ 13 kpc), even though these authors find little overall 
evidence for increased AGN activity in galaxy mergers. 

We caution that the correlation between AGN luminosity 
and projected separation and that of [Oiii] EW between the 
interacting components in an AGN pair support the hypoth- 
esis that galaxy tidal interactions enhance AGN activity, but 
do not directly test whether all AGNs are triggered in or as- 
sociated with galaxy interactions. Galaxy interactions may 
not be a necessary condition to trigger AGNs - secular pro- 
cesses may suffice at least for low-luminosity AGNs; nor is 
it a sufficient condition - there are internal conditions that 
need to be fulfilled such as the presence of ample quantities of 
gas. Whether galaxy interactions enhance AGN activity, and 
whether it is a sufficient condition to trigger AGNs are two 
separate issues. 

4.1.3. Is There A Time Delay between the Tidally-induced AGN and 
Star Formation Activity? 

To reconcile the elusive link between AGN activity and 
galaxy interactions found by some authors with the defini- 
tive connection between star formation and interactions, it 
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FIG. 12. — The fraction of galaxy pairs hosting AGNs (either single or 
double) among all galaxy pairs /agn (on a logarithmic scale) as a function of 
separation, /agn shows no significant dependence on r p from ~ 100 fr^, kpc 
to > 5 Kj l kpc for all AG Ns (shown in black), consistent with the clustering 
results of Li et al. (2006). On the other hand, for the stronger AGNs (e.g., 
L[Oin] > 10 40 erg s corresponding to LboI > 10 43 5 erg s~'), /agn rises 
prominently with decreasing r p for r p < 30 Mzk kpc (shown in blue). The 
fraction of galaxy pairs with both members hosting stronger AGNs /agn-agn 
also rises with decreasing r p for r p < 30 fr^ kpc (shown in red). Error bars 
indicate Poisson uncertainties. 

has been suggested that there is a considerable time delay 
(^0.1-1 Gyr) between th e interaction-induced star formation 
and AGN activity ( e.g.. ISanders et al.) 119881: iLi et al.ll2008bt 
lEllison et all 120081: fSchawins kTet allboiOl) For our AGN 



pairs, the values of r p below which the enhanced AGN and 
recent star-formation activity becomes significant are simi- 
lar, < 10-30 h^Q kpc. This suggests that there is no signif- 
icant time delay between the onset of interaction-enhanced 
AGN and star formation activity. However, as our sample 
does not probe the peak epoch of interaction-induced activ- 
ity, which is likely at later stages with much smaller nuclear 
separations than is probed by our sample, we cannot rule out 
a considerable time delay between the peak or overall aver- 
age epochs of tidally induced star formation and AGN activ- 
ity. A statistical sample of AGN pairs with smaller separa- 
tions is needed to better understand the tidal effects close to 
the peak epoch of the induced ac tivity. A complementary ap- 
proach to identifying AGN pairs dLiu et al.ll2010at IShen et ail 
1201 lat iMcGurk et al1l201 \i iFu et al.ll201 ll) based on the se- 
lectio n of AGNs with double-peaked narrow emission lines 
(e.gjHeckman et alj!981l:IZhou et al.l2 004; Com erford et al.l 
I2009t IWang et al.l 12009b iLiu et al.ll2010bt ISmith et al.ll2010h 
combined with follow-up spatially resolved near-IR (NIR) 
imaging and optical/NIR spectroscopy is more sensitive to 
smaller separation pairs (< a few kpc). 

4.2. Comparison with Simulations of Galaxy Mergers 

Our results are in broad agreement with the picture in 
which galaxy-galaxy interactions induce substantial gas in- 
flows which in turn promote star formation and AGN 
activity (e.g.. lNoguchi& Ishibashil 119861: iHernauistl 119891: 
iMihos & Hernquistl Il996t iPerez et al.l 120061) . although our 
sample only probes the low-luminosity regime. Current sim- 



ulations cannot resolve the scales relevant for star formation, 
AGN feeding, and stellar feedback, and thus employ simpli- 
fied recipes and assumptions to model star formation and the 
radiated AGN luminosity using galaxy properties on larger 
scales. Nevertheless, simulations suggest that, depending on 
progenitor properties and orbital parameters, the enhance- 
ment of star formation activity due to tidal effects may vary 
from levels as low as about twice the isolated case to values 
typical of starburst galaxies (20-60 times the isolated case) 
and the peak occurs during the interpenetrating phases (e.g., 
iDi Matteo et al.ll2007l) . The enhancement of recent star for- 
mation activity we observe (0.7-0.9±0.2 dex from r p > 100 
to ~ 5 h^Q kpc) is at the upper end of simulation predictions 
at the first close pas sage (1 to 4-5 times the isolated case; 
IDi Matteo et al.ll2007h . The enhancement in AGN luminos- 
ity we observe (~ 1 .0 dex in L[oni] from r p > 100 to ~ 5 h^ 
kpc; Figure|4]i is comparable to the prediction from numerical 
simulations (e.g., ~ 1 dex higher in SMBH acc retion rate just 
after th e first passage than the isolated case; ISpringel et al.l 
l2005bh . 

Numerical simulations of disk galaxies accreting dwarf 
satellites suggest that the primary disk galaxy develops a 
strong two-armed spiral pattern due to tidal perturbation from 
the infalling satellite, which drives sub stantial gas into the 
cente r and may fuel starburst ac tivity (M ihos & Hernquistl 
19941: iHernquist & Mihosl 1 19951) . On the other hand, 



Cox et al. I (120081) find very little, if any, induced star forma- 



tion for disk mergers with mass ratios larger than 5 : 1, es- 
pecially for galaxies with significant bulge components. The 
host galaxies in our samples of AGN pairs typically contain 
substantial bulge components, so our detection of enhanced 
recent star formation activity in AGN pairs with stellar mass 
ra tios M^i/M^ > 3 is in mild conflict with the prediction 
of ICox et alj d2008l) . This discrepancy could be due to differ- 
ences in the host-galaxy properties such as mass and struc- 
ture and interaction orbital properties. Our results are based 
on a special subset of all galaxy pairs, i.e., those in which 
both galaxies are AGNs. The simulations are not examing 
this subset exclusively, so that it is not a fair comparison. Al- 
ternatively it could point to the importance of physical pro- 
cesses neglected in simula tions (e.g., recycled gas from stel- 
lar winds and supernovae: ICiotti & Ostr iker 2007). Simula- 
tions with more tailored host-galaxy and interaction parame- 
ters may provide useful clues to discriminate between these 
possibilities. 

Simulations of disk galaxies with comparable masses sug- 
gest that the presence of dense central bulges delays the onset 
of gravitational instability, and gaseous inflows are weaker 
and o ccur earlier in bulgeless galaxies (Mi hos & Her nquist 
1996). Our result, that the star-formation enhancement is 
larger and becomes significant earlier (i.e., at larger r p ) in 
late-type than in early-type hosts of AGN pairs, broadly sup- 
ports this hypothesis. The observed difference is small, how- 
ever, probably due to the limited range of interaction phases 
spanned in our sample. Again, a statistical sample of AGN 
pairs with smaller separations is needed to enlarge the dy- 
namic range probed in separation. 

5. SUMMARY AND CONCLUSIONS 

We have studied recent star formation and BH accretion 
properties in the host galaxies of 1286 AGN pairs at z ~ 0.08 
and the subset of 256 pairs with unambiguous morphological 
tidal features, selected from the SDSS DR7 as described in 
Paper I. We have examined the effects of galaxy tidal interac- 
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tions on AGN pairs calibrated against control samples of iso- 
lated AGNs matched in redshift and stellar mass. We have in- 
vestigated the correlations between recent star formation and 
AGN activity and interaction parameters, the dependence of 
these correlations on host-galaxy properties, the correlation of 
star formation and AGN activity between the interacting com- 
ponents in each AGN pair, and the correlation between recent 
star formation and AGN activity in AGN pairs as compared 
to that observed in ordinary AGNs. Our main findings are 
summarized as follows. 

1 . We have found that the strengths of the 4000 A break 
and HS absorption in the host galaxies depend on the 
projected separation r p in AGN pairs with tidal features. 
They also depend weakly on LOS velocity offset Av. 
AGN pairs with smaller r p (Av) have smaller D„(4000) 
and larger H<5a, which indicate younger mean stellar 
ages and a higher fraction of post-starburst populations 
in the past 0.1 Gyr. Assuming the SSFR-D„(4000) cal- 
ibration of Brinchmannet al. (2004), the inferred me- 
dian enhancement of specific star formation rates is 
-0.7-0.9 ±0.2 dex from r p > 100 to - 5 h^ kpc. The 
enhancement of recent star formation becomes signifi- 
cant (> 3cr above the isolated case) at r p < 10-30 h^ 
kpc. 

2. We have found that the [Om] A5007 emission-line lu- 
minosity L[oni] and the Eddington ratio (using L[om] 
to infer bolometric luminosity and adopting stellar ve- 
locity dispersion as an indic ator for black hole mass 
assuming the calibration of iTremaine et al.l 120021) de- 
pend on r p in AGN pairs with tidal features, although 
the significance is smaller and with the trend being 
less prominent than that seen for recent star forma- 
tion. A weak dependence on Av is also detected. AGN 
pairs with smaller r p (Av) on average have larger Lpiii] 
and higher Eddington ratios, and the median enhance- 
ment of L[oni] (Eddington ratio) is — 0.7 ±0.1 dex 
(~ 0.5 ± 0.1 dex) from r p > 100 to - 5 h^ kpc. The 
enhancement of AGN activity becomes significant at 
r p < 10-30 h^Q kpc. This moderate overall increase 
we have detected in AGN luminosity/Eddington ratio 
in close AGN pairs is a result of tidally enhanced BH 
accretion events on top of a background of weak AGN 
events driven by secular processes not associated with 
interactions. 

3. We have compared how the effects of tidal interactions 
on AGN pairs vary when we are restricted to pairs with 
tidal features or not. We have found that the depen- 
dence of recent star-formation/AGN activity on r p is 
present in AGN pairs in both cases. On the other hand, 
the difference between the pair and control AGN sam- 
ples in their distributions of the examined BH accretion 
and host star formation properties is much smaller than 
that observed between the tidal and control samples. 
The discrepancy may be explained by contamination 
from closely separated pairs that are not tidally inter- 
acting; In addition, some AGN pairs without tidal fea- 
tures may indeed be interacting but are less affected by 
tidal effects, such as mergers involving more spheroidal 
components or those on retrograde orbits. Neverthe- 
less, our results suggest that the effects of galaxy tidal 
interactions will be underestimated using overall prop- 



erties of galaxy pairs selected solely based on pro- 
jected separation and velocity offsets, in particular if the 
adopted pair separation threshold is substantially larger 
than 10-30 h^ kpc. 

4. We have explored how host-galaxy mass ratio and 
structure regulate the effects of tidal interactions in 
AGN pairs. We have found that, at 5 < r p < 100 h^,' 
kpc, the relative enhancement level of star-formation 
activity in interacting AGN pairs with stellar mass ra- 
tios M*.i/M*.2 > 3 is comparable to those with com- 
parable stellar masses (M^ j/M*^ < 3) in both com- 
ponents. While the enhancement of AGN activity is 
observed in AGN pairs of stellar mass ratios as high 
as 10, the relative enhancement on average is more 
prominent in the secondary than in the primary AGN 
(e.g., 1.2 ± 0.3 dex compared to 0.5 ± 0.4 dex in the 
median value of L[onrj> or 0.7 ±0.3 dex compared to 
0.6 ± 0.4 dex in the median value of Eddington ratio). 
We have detected no significant difference in the tidally 
enhanced activity when we divide the host galaxies ac- 
cording to their r-band concentration index as a proxy 
for morphology, but there is a marginally larger en- 
hancement for late-type hosts. In addition, there is ten- 
tative evidence that the threshold r p value below which 
the recent star-formation enhancement becomes signif- 
icant is larger in late- than in early-type hosts. These re- 
sults may re flect the stabilization effect o f bulges as pre- 
dicted by the lMihos & Hernquistl (H996) galaxy merger 
simulations. 

5. We have observed a statistically significant correla- 
tion of recent star formation (indicated by H<5a and 
D„(4000) ) between the interacting components in an 
AGN pair (albeit with a large scatter), similar to the 
color correlation known in inactive galaxy pairs. Simi- 
larly, we have detected a statistically significant correla- 
tion of AGN activity as indicated by [Om] A5007 rest- 
frame equivalent width between the two components in 
an AGN pair, although the scatter is significant. 

6. We have found a correlation between recent star forma- 
tion and AGN luminosity in AGN pairs similar to that 
observed in control samples of ordinary AGNs. While 
AGN pairs follow the same scaling relation of control 
AGNs, their host galaxies exhibit younger mean stellar 
ages, higher starburst fractions, and host more powerful 
AGNs. 

In this paper, we have examined the properties of AGN 
pairs and compared them with control samples of single 
AGNs to characterize the effects of tidal interactions. The 
effects that we have found are based on a special population 
of galaxy pairs in which the two central SMBHs are both ac- 
tive at the same time. More generally, galaxy pairs often host 
single or no AGNs. Are galaxy pairs hosting double AGNs 
intrinsically different from those hosting single or no AGNs 
(e.g., in terms of host galaxy, BH properties, and/or local or 
large-scale environments), or is the double-AGN phase uni- 
versal to all galaxy tidal interactions? Is the triggering of both 
AGNs preferential to certain merger configurations, or are the 
two BHs just being turned on and off stochastically? In a sub- 
sequent paper, we will compare AGN pairs with galaxy pairs 
hosting single or no AGNs to answer these questions. 
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The analysis of z ~ 0.1 AGN pairs as a function of sep- 
aration or velocity offset as a proxy for merger stage pre- 
sented here also needs to be extended to smaller pair sepa- 
rations (i.e., more advanced stages) and to higher redshifts, 
where the effects of mergers and close interactions are pre- 
dicted to be even more dramatic. Recently there has been an 
increasing number of AGN pairs with r p < a few kpc iden- 
tified based on the selecti on of AGNs with double-peaked 
narrow emission lines (e.g., [H eckman et al. 1981: [Zhou et al. | 
2004 IComerford e t al ] 12 009: Wa ng et al.n2009l; iLiu etaLl 
2010b; ISmith et all 1201 Oh combined with follow-up spa- 
tially reso lved near-IR (NIR) imaging and op t ical/NIR spec- 
troscopy (ILiu et al] 1201 Oat iShen et all 120711 IMcGurk et al.1 
1201 U iFu et all 1201 ll) to discriminate between double AGNs 
and single AGNs with complex gas kinematics. This ap- 
proach in principle can identify AGN pairs with separations 
as small as ~ 100 pc, i.e., down to the limit set by the intrin- 
sic size of narrow line regions around si ngle AGNs. At z > 1 , 
i.e., close to the peak AGN epoch (e.g.. iMarconi et aLll2004t 
Hasin ger et al] 12005). there have b een studies on the local 
environments of single AGNs (e.g.. | G eorgakakis et al .1120081: 
Montero-Dorta et al j|2009l: ICoil et a.1.1 12009b ISilverman et al] 
2009] |2011|) in the DEEP2 survey dDavis et al.1 120031) or the 
zCOSMOS survey dLillv et al.ll2007blScoville et al.ll2007l) us- 
ing X-ray identification of AGNs based on Chandra obser- 
yations in the A EGIS dDavis et all 120071) or C-COSMOS 
dElvis et all 120091) field, although a statistical study of AGN 
pairs similar to the analysis present in this paper is still ham- 
pered by small galaxy sample size combined with the diffi- 
culty of resolving close pairs and of identifying tidal distur- 
bance features. Future multi-object NIR sp ectrograph such 
as the MOSFIRE/Keck (McLean et alll2008l) combined with 
high resolution and deep imaging capability such as offered 
by the HST or ground-based adaptive optics facilities will en- 
able observations of galaxy samples with sufficient statistical 
power, resolution, and sensitivity to perform a similar study 
of AGN pairs at z > 1 . 
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